Direct observation of sub-binomial light 
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Nonclassical states of light are necessary resources for quantum technologies such as cryptography, compu- 
tation and the definition of metrological standards. Observing signatures of nonclassicality generally requires 
inferring either the photon number distribution or a quasi-probability distribution indirectly from a set of mea- 
surements. Here, we report an experiment in which the nonclassical character of families of quantum states is 
assessed by direct inspection of the outcomes from a multiplexed photon counter. This scheme does not register 
the actual photon number distribution; the statistics of the detector clicks alone serve as a witness of nonclas- 
sicality, as proposed by Sperling et al in Phys. Rev. Lett. 109, 093601 (2012). Our work paves a way for the 
practical characterisation of increasingly sophisticated states and detectors. 



States of the electromagnetic field with no analogue in a 
classical theory of electromagnetism provide a means to char- 
acterise quantum coherence and are a valuable resource for 
quantum-enhanced technologies. Determining what consti- 
tutes a nonclassical state not a trivial matter. Many crite- 
ria to establish where the quantum-classical boundary lies 
have been proposed. Typically, one considers pathologi- 
cal behaviours of distributions in phase space: one looks 
for negative values of the Wigner quasi-probability function 
on the quadrature space or values more singular than a 
delta function of the Glauber-Sudarshan P distribution on the 
phase-amplitude space [T-^, with the latter definition cover- 
ing a larger class of states. When detecting photons directly, 
the Mandel parameter Qm is often used 0. This neatly cap- 
tures any sub-Poissonian behaviour in the photon statistics as 
a witness of nonclassicality. 

Experimentally, these criteria have been widely adopted, 
and inferred either by full reconstruction of the distribu- 
tion L6rrl0iL or by direct measurement [TT] [121. A differ- 
ent approach consists in examining the photon number dis- 
tribution in a counting experiment 1 5] [TsUTtI l25ll . While ac- 
tual photon number resolving detectors are becoming more 
frequently available |[T9ti22ll . simpler solutions are still ap- 
pealing as quantum networks become more complex. Mul- 
tiplexed detectors are a commonplace choice for accessing 
higher order Fock states by binning the light to several on- 
off detectors |^[24l. As such, these detectors do not out- 
put the photon number distribution. It is recognised that the 
Mandel parameter applied to the click distribution would give 
false indications of nonclassicality even for a large number 
of detection bins 1.27 J . A linear inversion technique is then 
required |[IIl[l8||26l|33 in order to obtain the actual distribu- 
tion from the distribution of the clicks. However, this becomes 
more sensitive to noise and less tractable with increasing sys- 
tem size. 

In this paper, we directly observe a signature of nonclassi- 
cality in the output of a time-multiplexed detector, following 
the proposal in Ref. |28 1 which identifies a sub-binomial be- 
haviour in the click statistics. We experimentally measure this 
sub-binomial character for three families of quantum states 



and compare it with the standard Mandel parameter L5J, as 
well as a naive analysis of the data. Our results show that the 
sub-binomial behaviour of click statistics is a reliable indica- 
tor of nonclassicality of single-mode light fields. Moreover, 
this test performs just as well at inferring the true Mandel pa- 
rameter from the photon distribution. As quantum informa- 
tion systems evolve in sophistication, characterisation via full 
tomography becomes prohibitively expensive, and our results 




FIG. 1: (Color online) a: General scheme of a multiplexed detec- 
tor. Access to more detail of the photon distribution is obtained by 
dividing the input beam onto several modes, which are then mea- 
sured with on/off detectors, such as APDs. The photon number dis- 
tribution {pn} determines the statistics of the counts {c^}, i.e. the 
probabilities of observing events with i clicks regardless the bin in 
which they occur, b: Spatially multiplexed detector: the input beam 
is split by a set of beam splitters onto a collection of spatial modes, 
c: Time-multiplexed detector: the input beam is coupled in a single- 
mode fibre and split on a first 50:50 beam splitter. The two arms 
are then recombined to form an unbalanced Mach-Zehnder interfer- 
ometer with a relative delay At=50ns. This produces four distinct 
modes: two time bins on two distinct spatial modes. This operation 
is repeated in a second unbalanced interferometer with double rela- 
tive delay 2At, in order to produce two more time bins. This is the 
detector we adopted in the present experiments. 
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pave the way for a simple characterisation of quantum states 
and detectors. 

The general scheme of a multiplexed detector is illustrated 
in Fig. la. These simple devices permit access to more detail 
about the number distribution of a light field through sequen- 
tial division of an input beam on a series of beam splitters, 
with detection of each output mode by avalanche photodiodes 
(APDs). Different architectures have been realised by using 
either spatial (Fig. lb) | 35 | or time-bin (Fig. Ic) | 36| multi- 
plexing. In this way, APDs with no photon number resolving 
capabilities can still deliver information about an input state 
with n>l photons. However, this is non-deterministic, occur- 
ring only when the photons split into TV separated spatial or 
temporal bins. This means that the distribution of the pho- 
tons {pn}, and the one of the clicks from the APDs {c^} are 
different, and their difference only scales with the number of 
bins as 1/N |27|. Some characterisation of the device is re- 
quired to find the matrix relating the two distributions. Since 
the splitting operates a linear transformation, the matrix can be 
inverted in order to obtain the values of {pn} fT\},^24^,^29\. In 
our realisation, we adopt a time-multiplexed detector (TMD) 
which can detect up to 7V=8 events by splitting the input 
into two spatial and four temporal modes, whose separation 
is larger than the typical dead time of an avalanche photodi- 
ode l[T8ll23l [2411261. 

The question arises whether there exists a way of observing 
nonclassicality directly in the measured quantities {ci}. A 
solution can be found by inspecting the Mandel parameter: 



n 



-1. 



(1) 



The negativity of this parameter establishes a suffi- 
cient criterion for nonclassicality: this implies that 
the mean n= Xl^o is larger than the variance 
A^n=^^Q(n — n) Pn, i.e. if the statistics are sub- 
Poissonian. Coherent states, which benchmark classical light, 
display Poissonian number statistics. By contrast, Fock states 
display sub-Poissonian statistics, hence the corresponding 
Mandel parameter takes negative values. When using TMDs, 
a Poissonian distribution of the photons {pn} results in a bi- 
nomial distribution of the counts {ci} [27, .28]. Based on this, 
Sperling, Vogel and Agarwal (SVA) have developed a differ- 
ent nonclassicality criterion applying to the measured click 
statistics: 
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where c = Y^f^QiCi, and A^c=^^Q(i — c)^q. This ex- 
tends Mandel 's formula ([T]) to the case of a binomial dis- 
tribution with mean x, the variance of which is given by 
x{l — x). Consequently, nonclassicality can be directly in- 
ferred from the sub-binomial behaviour of the observed statis- 
tics, i.e. Qb < 0. 

Here, we experimentally measure both Qb and Qm for 
families of quantum states that may be tuned continually be- 
tween the classical and nonclassical regimes. To obtain such 
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FIG. 2: (Color online) a: Measurement of the nonclassicality of 
a single-mode state. Photons are produced by pulsed type-II co- 
linear parametric down-conversion in a KDP crystal driven by a 
doubled Ti:Sa pump laser (Ap=415nm, AA =50 nm, repetition 
rate 256 kHz). This produces pairs of in a single spatio-temporal 
mode |31 1. By heralding the presence of a photon in the V mode 
reflected by PBS i is heralded by a click of an APD on the H mode, 
and a trigger derived from the driving laser. A coherent state \a), 
derived from the main laser, is made to interfere with the single pho- 
ton on a variable beam splitter consisting of PBS2, a half- wave -plate 
HWP, and PBS3. The effective reflectivity of the device is then set 
by the angle of the HWP as = cos^ 0. The transmitted arm is 
delivered to detector TMDi which provides a second herald event, 
either no clicks A:=0 or a single click k=l. The conditional states on 
the reflected arm are then analysed by TMD2. b: Observing sub- 
binomial behaviours in the joint click statistics from a two-mode 
squeezed state generated from the same crystal at higher pump in- 
tensity. The two modes from the KDP are directly delivered onto 
two detectors TMDi and TMD2, and clicks are collected. 



states, we interfere a single photon with a coherent state 
on a beam splitter with variable reflectivity R, and consider 
the state on the transmitted arm conditioned on the presence 
of k photons on the reflected arm. For k=l, we produce 
non-Gaussian states by performing photon catalysis |26l |30l . 
These allow us to investigate the behaviour of in the ex- 
periment, as we tune the state from a single photon (R=l) to 
a coherent state (R=0). 

The experimental scheme is shown in Fig. [2]3: a single 
photon is produced by pulsed parametric down-conversion 
in a nonlinear KDP crystal. This interaction produces a 
pair of photons with orthogonal polarisations in a spectrally- 
factorable state |31 1: the horizontally (H) polarised photon is 
detected by an APD and serves as a herald event. The verti- 
cally polarised (V) photon is spatially overlapped with a co- 
herent state I a), i^-polarised, on a polarising beam splitter 
(PBS). A variable beam splitter, comprising a half- wave-plate 
and a second PBS, is then used to interfere the two beams. 
We monitor one output port ,heralding on the measurement 
of either k=0, or k=l photons by a first TMD. The condi- 
tional state on the second port constitutes our quantum signal, 
and the complete click statistics from a second TMD {ci} are 
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collected. 

We perform three different statistical analyses on the col- 
lected data, whose results are reported in Fig. [3^ for k=0, and 
in Fig. (Sj) for k=l. In the first analysis, we assess the sub- 
Poissonian character of the click statistics using a Mandel-like 
parameter (red squares) defined as: 
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The blue dots report the experimental result for the genuine 
Mandel parameter Eq. ([T]), after the photon distribution has 
been obtained by inversion 1 37 1. Finally, the experimental val- 
ues for the binomial parameter Eq. ^ are shown as the black 
squares. 
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FIG. 3: (Color online) Test of nonclassicality for single-mode states 
heralded by (a) /c=0, and (b) k=l. In both figures, the black boxes 
show the the Mandel parameter Qm and the blue points show the 
S VA parameter Q b • The values of the Mandel-like parameter calcu- 
lated from the clicks Qf are shown as red points. The limit ^ 
correspond to coherent states for which one expects no nonclassical- 
ity, while around R ^ 1, the states approximate displaced single- 
photons |26|. All experimental points are obtained by 600 s of data 
collection. Errors are evaluated by a Monte Carlo technique assum- 
ing Poissonian noise. 

In general, the Mandel-like parameter Qp cannot be used 
to detect nonclassicality, as it can take negative values even 
for the coherent states obtained when i?~0 L27J . as shown in 
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FIG. 4: (Color online) a. Joint statistics for detection with TMD 
on two-mode squeezed states, due to limited efficiency (^ 7%), the 
most frequent event is when no clicks are registered at either side. All 
the events in coincidence with a trigger signal from the main laser are 
included; b. Measured Q b parameter for the state detected by TMD 2 
conditioned on different click events on TMDi c. Statistical analysis 
of the Q B parameter when conditioning on two clicks from TMD2 
(left), and TMDi (right). The histograms are obtained based on 10^ 
Monte Carlo simulations assuming Poissonian noise. 



Fig. [3}). Further, its behaviour depends on the particular class 
of states under investigation: for the case in Fig. [3^ the sepa- 
ration between this parameter and the Mandel parameter can 
be significant, while they are closer in Fig.|3|3. In both cases, 
a systematic failure of Qf can be seen at low reflectivities. 
This issue is resolved by employing the binomial parameter 
Qb, which performs as well as the Mandel parameter Qm for 
our eight-bin detector, and with similar uncertainties. In our 
experiment, we have investigated an intensity regime in which 
the inversion is not affected by the high-order time bins, in or- 
der to have a reliable quantity of comparison for the SVA pa- 
rameter of our states. However, this might not be feasible in 
the general case, as uncertainties in these terms will propagate 
badly in the reconstructed {pn}- Thus the direct approach of 
SVA is useful in avoiding this source of error, as well instabil- 
ities arising from noise. 

We also use the SVA parameter to characterise the output 
from the down-converter at higher gain, producing a two- 
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mode squeezed state of the form |^)(x J2n>o ^^\^)h\i^)v on 
the two orthogonally polarised output modes. If we were to 
use a photon-number resolving detector, the observation of n 
photons on the H mode would imply the presence of n pho- 
tons, hence a sub-Poissonian state, on the V mode 1321 1331 . 
In our experiment, as shown in Fig. 2b, we observe a similar 
behaviour in terms of sub-binomial light: a click from a TMD 
projects the other mode on to a sub-binomial state. The joint 
statistics from two TMDs, each measuring a different polari- 
sation mode is shown in Fig.|4^: from this we have performed 
the S VA test on the states obtained by conditioning on a given 
number of clicks on one TMD or the other. 

When all the events of TMD2 are considered regardless 
of the measurement result at TMDi, or vice versa, we ob- 
serve classical behaviour, as expected for a thermal state 
(Fig.|4]3): we obtain the values Qb = (9.3 ± 0.6) • 10"^ 
for the state on mode 1, and = (10.9 ± 0.6) • 10"^ 
for the state on mode 2. Discrepancies with respect to the 
expected symmetry can be attributed to the different detec- 
tion efficiencies which have been measured. When only the 
events causing no clicks on TMDi are registered, we still ob- 
serve a similar behaviour: the reduced efficiency of our detec- 
tor, in fact, makes this event similar to tracing out the other 
mode. The measured values are = (10.2 ± 0.6) • 10~^ 
(mode 1) and Qb = (8.5 ± 0.6) • 10"^ (mode 2). When 
single clicks are considered as heralds, the state becomes sub- 
binomial, and so it happens independently of the mode we 
use as the trigger: Qb = (-3.84 ± 0.33) • 10"^ (mode 1), 
and Qb = (-4.49 ± 0.30) • 10"^ (mode 2). Finally, when 
triggering on two clicks, the state does present some evidence 
for nonclassicality, Qb = (—3.3 ± 3.9) • 10~^ (mode 1) and 
Qb = (-8.3 ± 3.5) • 10"^ (mode 2), although the low counts 
allow us to be confident only for mode 2. A statistical analy- 
sis conducted by simulating experiments using a Monte Carlo 
technique is reported in Fig. [4]:: while the most likely values 
for both modes do lie well below zero, the tails of the distri- 
bution make the value of Qb for mode 1 still compatible with 
positive values. 

Our experimental results, in the light of the retrodictive ap- 
proach in 1 34 1, also provides preliminary evidence for a test 
of nonclassicality of the TMDs themselves. In the retrodictive 
framework, the nonclassicality of a detector D can be ascer- 
tained by observing the nonclassicality of one mode of |?/^) 
when the other mode is projected by the action of D in the 
limit of high squeezing A^l. Although we are far from this 
regime, the observed negativity of Qb resulting from condi- 
tioning one of the arms on a TMD points towards the nonclas- 
sicality of the device. 

The direct experimental observation of the nonclassicality 
of light fields is a vital aspect of foundational studies and tech- 
nological applications of quantum mechanics. In this work 
we have shown how this is made possible by the use of time- 
multiplexed detectors, and by directly inspecting the resulting 
click statistics. This represents a more elegant and unambigu- 
ous procedure for confirming quantumness than their indirect 
verification by deconvolving the action of the detector. 
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